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Medium modification of parton fragmentation functions induced by multiple scattering and gluon 
bremsstrahlung is shown to describe the recent HERMES data in deeply inelastic scattering (DIS) 
very well, providing the first evidence of ^4 2 / 3 -dependence of the modification. The energy loss is 
found to be (dE/dL) ~ 0.5 GeV/fm for a 10-GeV quark in a Au nucleus. Including the effect of 
expansion, analysis of the ty° spectra in central Au + Au collisions at y^i = 130 GeV yields an 
averaged energy loss equivalent to (dE/dL) ~ 7.3 GeV/fm in a static medium. Predictions for 
central Au + Au collisions at -y/s = 200 GeV are also given. 
PACS numbers: 12.38. Mh, 24.85. +p; 13.60.-r, 25.75.-q 



Energetic partons produced via hard processes pro- 
vide an excellent tool that enables tomographic studies 
of both hot dense and cold nuclear matter. By measur- 
ing the attenuation [1-5] of these partons as they prop- 
agate through the medium, one would be able to study 
the properties such as the geometry [6,7] and the gluon 
density of the medium. The attenuation will suppress 
the final leading hadron distribution giving rise to modi- 
fied parton fragmentation functions [8]. Such a modified 
fragmentation function inside a nucleus has been derived 
recently [9] in a pQCD approach with a systematic ex- 
pansion of higher-twist corrections to the fragmentation 
processes. In this letter, we will compare the predicted 
nuclear modification to the recent HERMES experimen- 
tal data [10] and extract the effective parton energy loss. 
We will further extend the study to the case of a hot 
QCD medium including the dynamics of expansion. We 
will then analyze within this framework the 7r° spectra 
as measured by the PHENIX experiment [11] in central 
Au + Au collisions at i/s = 130 GeV, which have shown 
significant suppression at large transverse momentum. 
The extracted effective parton energy loss will be com- 
pared with that in a cold nucleus, and discussions will 
be given about the implications of the PHENIX data on 
the gluon density in the early stage of central Au + Au 
collisions at the RHIC energy. We will also provide pre- 
dictions for 7T° spectra in central Au + Au collisions at 
y/s = 200 GeV. 



In deeply inelastic scatterings (DIS), quark fragmen- 
tation functions are factorizable in leading twist from 
the parton distribution functions and photon-quark (7*9) 
scattering cross section. In a nucleus target, the quark 
struck by the virtual photon suffers multiple scattering 
and induced bremsstrahlung before hadronization. Ex- 
tending the generalized factorization [12] to the semi- 
inclusive process, e(Li) + A(p) — > e(L 2 ) + h(£h) + X, 
one can define a modified fragmentation function as [9], 
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where f q {x, Q 2 ) is the quark distribution function in the 

nucleus, s = (p + L\) 2 and p — [p + , 0, OjJ is the momen- 
tum per nucleon inside the nucleus. The momentum of 
the virtual photon 7* is q = [—Q 2 /2q~,q~,0±] and the 
momentum fraction carried by the hadron is Zh = £Z I q~ ■ 

The hard part of the 7*5 scattering HfrJ is the same as 
in ep scattering [9] and L M „ = \ Tr(^i7 M ^Iv)- 

Including the leading twist-4 contributions from double 
scattering processes, the modified effective quark frag- 
mentation function can be obtained as [9] 

D q ^ h {z hl Q 2 ) = D q ^ h {z h ,Q 2 ) + AD q ^ h (z h ,Q 2 )-, (2) 
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[A'J q ^ q g(z, X, XL,£T)Dq^h(Zh/z) + Aj q ^ gq (z, X,XL,^T) D 9^h(Zh/z)] , (3) 



A-i q ^ q g(z,x,x L ,t^) = 

where, Ca = 3, N c = 3, A r y q ^ gq (z,x,xi J ,£ 2 ^) = Aj q ^ qg (l — z,x,xl,£t) are the modified splitting functions, xl = 
l\ jlp^ q~ z(A — z), and D a ^h{zh, Q 2 ) are the normal twist-2 parton fragmentation functions in vacuum. The (5-function 
part in the modified splitting function is from the virtual corrections, with AT^(i,4) defined as 
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AT A (xJ 2 T ) = f dz-^—[2T A (x,x L )\ z=1 -(l + z 2 )T A (x,x L )]. (5) 

Such virtual or absorptive corrections are important to ensure the infrared safety of the modified fragmentation 
function and the unitarity of the gluon radiation processes. The quark-gluon correlation function 

T A (x,x L ) = j ^ dy^ dy^ e i(x+XL)p+v ~ +iXTP+{y ^ - v z ] {1 - e r lXLV+ ^){\ - e ~ ix ^ + iy~ ~vT)) 

\ {A\4, q {Q) 1 +F+{y^)F + °(y^)iP q {y-)\A)9{-y^)6(y-- y ^) (6) 



contains essentially four independent twist-4 parton ma- 
trix elements in a nucleus [xt — (k^) /2p + q~ z(l — 
z)]. The dipole-like form-factor (1 - e - iXLP+v ^){\ - 
e -ix L p+(y -y 1 )^ ar i ses f rom the interference between the 
final state radiation of the 7*5 scattering and the gluon 
bremsstrahlung induced by the secondary quark-gluon 
scattering. By generalizing the factorization assumption 
[12] to these twist-four parton matrices, we have 

T q A g (x,x L ) » C(Q 2 )m N R A f A (x)(l - e"^), (7) 

with a Gaussian nuclear distribution p(r) <~ 
exp{-r 2 /2R 2 A ), R A = 1.12A 1 / 3 fm. Here, x A = 
I/ttinRa, and is the nucleon mass. We should 
emphasize that the parameter C(Q 2 ) should in principle 
depend on the renormalization scale Q 2 among other 
kinetic variables as shown recently in a detailed analysis 
of the twist-four nuclear matrix elements [13]. Therefore, 
it can take different values in different processes. 
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FIG. 1. Predicted nuclear modification of jet fragmenta- 
tion function is compared to the HERMES data [10] on ratios 
of hadron distributions between A and D targets in DIS. 

Since the two interference terms in the dipole-like 
form- factor involve transferring momentum xlP + be- 
tween different nuclcons inside a nucleus, they should 
be suppressed for large nuclear size or large momen- 
tum fraction xl- Notice that r/ = 1/xlp + is the 
gluon's formation time. Thus, xl/x a — La/tj, with 
La = Ra^n /p + being the nuclear size in the infinite mo- 
mentum frame. The effective parton correlation and the 



induced gluon emission vanishes when the formation time 
is much larger than the nuclear size, xl/xa <C 1, because 
of the Landau-Pomeranchuck-Migdal (LPM) interference 
effect. Therefore, the LPM interference restricts the ra- 
diated gluon to have a minimum transverse momentum 
£?r ~ Q 2 /m N R A ~ Q 2 /A 1 / 3 . The nuclear corrections to 
the fragmentation function due to double parton scatter- 
ing will then be in the order of a. s A x l 3 jl 2 T ~ a s A 2/3 /Q 2 , 
which depends quadratically on the nuclear size. For 
large values of A and Q 2 , these corrections are leading; 
yet the requirement £ 2 r <C Q 2 for the logarithmic approx- 
imation in deriving the modified fragmentation function 
is still valid. 
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FIG. 2. Energy dependence of the nuclear modification 
compared with the HERMES data [10]. 



With the assumption of the factorized form of the 
twist-4 nuclear parton matrices, there is only one free 
parameter C(Q 2 ) which represents quark-gluon correla- 
tion strength inside nuclei. Once it is fixed, one can pre- 
dict the z, energy and nuclear dependence of the medium 
modification of the fragmentation function. Shown in 
Figs. 1 and 2 arc the calculated nuclear modification 
factor of the fragmentation functions for 1A N and s4 Kr 
targets as compared to the recent HERMES data [10]. 
There are strong correlations among values of Q 2 , v and 
z in the HERMES data which are also taken in account 
in our calculation. The predicted shape of the z- and 
independence agrees well with the experimental data. 
A remarkable feature of the prediction is the quadratic 
A 2 / 3 nuclear size dependence, which is verified for the 



first time by an experiment. This quadratic dependence 
comes from the combination of the QCD radiation spec- 
trum and the modification of the available phase space 
in It or xl due to the LPM interferences. 

The observed attenuation of the leading hadrons could 
also be attributed phenomenologically to hadron absorp- 
tion inside the nucleus. This is, however, achieved only 
via some ad hoc assumption of the hadron formation time 
[10]. Considering the hadronization process as regenera- 
tion of gluon field within a spatial region of a hadon size 
r/j, the hadron formation time for a light quark will be 
r\ [14]. Taking w 1 fm, th is about 40 fm for 



t 



f 



v = 8 GeV, the lower limit of the HEREMS experiment. 
This is much larger than the size of the heaviest nuclei 
available. We therefore assume that the nuclear attenu- 
ation in this energy region is mainly caused by induced 
gluon radiation and multiple parton scattering. 

By fitting the overall suppression for one nuclear tar- 
get, we obtain the only parameter in our calculation, 
C(Q 2 ) = 0.0060 GeV 2 with a s (Q 2 ) = 0.33 at Q 2 w 3 
GeV 2 . This parameter is also related to nuclear broaden- 
ing of the transverse momentum of the Drell-Yan dilep- 
ton in pA collisions [15], (Ag 2 ) « C-Ka s /N C XA- With 
an experimental value of (Aq\) = 0.016A 1 / 3 GeV 2 [16] 
and a s (M 2 -) = 0.21 (< M 2 - 40 GeV 2 ), one finds 

<5(M 2 -) = 0.013 GeV 2 , which is about a factor 2 larger 
than the value obtained in our fit to the HERMES 
data. The value of C determined from nuclear broad- 
ening in photo-production of di-jet is even larger [12] at 
Q 2 = Ap\ rs 64 GeV 2 . Such a strong scale dependence of 
C(Q 2 ) is in line with one's expectation since it is related 
to gluon distribution xg(x, Q 2 ) at small x in nuclei [13]. 

If one defines theoretically the quark energy loss as 
that carried by the radiated gluons, then the averaged 
total fractional energy loss is [cf. Eq. (3)], 
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the averaged total energy loss as AE 

l2\„,2/n2\™..E>2 



v(Az g ) « 

C(Q 2 )a 2 {Q 2 )m N R\{C A /N c )Z\n{l/2x B ). With the de- 
termined value of C, (x B ) « 0.124 in the HERMES ex- 
periment [10] and the average distance (La) = Ra \/2 Jtt 
for the assumed Gaussian nuclear distribution, one gets 
the quark energy loss dE/dL w 0.5 GeV/fm inside a Au 
nucleus. 

To extend our study of modified fragmentation func- 
tions to jets in heavy-ion collisions, we assume (k T ) w p? 
(the Debye screening mass) and a gluon density profile 



p(y) = ( T o/ T )d(R A — y)po for a 1-dimensional expand- 
ing system. Since the initial jet production rate is in- 
dependent of the final gluon density which can be re- 
lated to the parton-gluon scattering cross section [2] [ 
a s XTG(xr) ~ p 2(J g\i one nas then 



a s T^ g {x B ,x L ) 
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[i 2 J dyo- g p(y)[l - cos(j//r/)], (10) 



where r/ = 2Ez{\ — z)/£ T is the gluon formation time. 
One can recover the form of energy loss in a thin plasma 
obtained in the opacity expansion approach [2,7], 
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Keeping only the dominant contribution and assuming 
(7 g w C a 2ira 2 /p 2 (C a — 1 for qg and 9/4 for gg scatter- 
ing), one obtains the averaged energy loss, 



,dE irC a C A a 3 rRA 
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Neglecting the logarithmic dependence on r, the aver- 
aged energy loss in a 1-dimensional expanding system 
can be expressed as 



x ~dL 
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dE 2t 
dL Ra 



(13) 



where dE /dL cx poRa is the energy loss in a static 
medium with the same gluon density p as in a 1-d ex- 
panding system at time To- Because of the expansion, 
the averaged energy loss (dE/dL) id is suppressed as com- 
pared to the static case and does not depend linearly on 
the system size. This could be one of the reasons why 
the effect of parton energy loss is found to be negligible 
in A A collisions at y/s = 17.3 GeV [20]. 

An effective model of modified fragmentation functions 
was proposed in Ref. [8] : 
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with Az to account for the fractional parton energy loss. 
This effective model is found to reproduce the pQCD re- 
sult from Eq. (3) very well, but only when Az is set to 
be Az w 0.6(z g ). Therefore the actual averaged parton 
energy loss should be AE/ E = l.QAz with Az extracted 
from the effective model. The factor 1.6 is mainly caused 
by unitarity correction effect in the pQCD calculation. 
A similar effect is also found in the opacity expansion 
approach [17]. 

The PHENIX experiment has reported [11] a strong 
suppression of high px hadrons in central Au + Au col- 
lisions at y/s = 130 GeV. To extract the parton energy 



loss, we compare the data with the calculated hadron 
Pt spectra in heavy-ion collisions using the above effec- 
tive model for medium modified jet fragmentation func- 
tions [18]. Shown in Fig. 3 are the nuclear modifica- 
tion factor Raa(pt) as the ratio of hadron spectra in 
AA (pA) and pp collisions normalized by the number 
of binary collisions [19]. Parton shadowing and nuclear 
broadening of the intrinsic fcy are also taken into account 
in the calculation which decrsibes pA data for energies 
up to y/s = 40 GeV [18]. The nuclear fc^-broadcning 
gives the so-called Cronin enhancement at large pt in 
pA collisions, where there is no parton energy loss in- 
duced by a hot medium. Fitting the PHENIX data yields 
(dE/dL) ld w 0.34 In Ej In 5 GcV/fm, including the fac- 
tor of 1.6 from the unitarity correction effect. We con- 
sider only 7T° data here, since at large px the charged 
hadrons are dominated by baryons, which could be influ- 
enced mainly by non-perturbative dynamics [7]. 
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FIG. 3. Calculated nuclear modification factor of 7r pr 
spectra for d + Au and central Au + Au collisions at y/s =130 
(solid) and 200 GeV (dashed) as compared to PHENIX 
data [11] (solid circles are measured spectra normalized by 
PHENIX parameterization of the pp spectra while open cir- 
cles are normalized by our calculated pp spectra). 

Taking into account the expansion, the averaged par- 
ton energy loss extracted from the PHENIX data would 
be equivalent to (dE/dL) = 0.34(.Ra/2to) ln.E/ln5 in 
a static system with the same gluon density as the ini- 
tial value of the expanding system at r . With Ra ~ 6 
fm and tq ~ 0.2 fm, this would give (dE/dL)o « 7.3 
GeV/fm for a 10-GeV parton, which is about 15 times of 
that in a cold Au nucleus as extracted from the HERMES 
data. Since the parton energy loss is directly proportional 
to gluon density of the medium, we can predict the ir° 
spectra at y/s — 200 GeV as given by the dashed lines in 
Fig. 3, assuming that the initial parton density in central 
Au + Au collisions at y/s = 200 GeV is about 10% higher 
than at 130 GeV. 



In summary, the nuclear modification of parton frag- 
mentation function predicted in a pQCD study describes 
well the HERMES experimental data. The extracted 
energy loss is dE/dL w 0.5 GeV/fm for a quark with 
E = 10 GeV in a Au nucleus. Analysis of the PHENIX 
data of 7T° spectra suppression in central Au + Au colli- 
sions yields an averaged parton energy loss in an expand- 
ing system that would be equivalent to (dE/dL) n w 7.3 
GeV/fm in a static medium. 
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